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ABSTRACT 

Although A^-body simulations of cosmic structure formation suggest that dark matter halos have 
density profiles shallower than isothermal at small radii and steeper at large radii, whether observed 
galaxy clusters follow this profile is still ambiguous. We use one such density profile, the asymmetric 
NFW profile, to model the mass distributions of 11 galaxy clusters with gravitational arcs observed 
by HST. We characterize the galaxy lenses in each cluster as NFW ellipsoids, each defined by an 
unknown scale convergence, scale radius, ellipticity, and position angle. For a given set of values of 
these parameters, we compute the arcs that would be produced by such a lens system. To define 
the goodness of fit to the observed arc system, we define a function encompassing the overlap 
between the observed and reproduced arcs as well as the agreement between the predicted arc sources 
and the observational constraints on the source system. We minimize this to find the values of 
the lens parameters that best reproduce the observed arc system in a given cluster. Here we report 
our best-fit lens parameters and corresponding mass estimates for each of the 11 lensing clusters. 
We find that cluster mass models based on lensing galaxies defined as NFW ellipsoids can accurately 
reproduce the observed arcs, and that the best-fit parameters to such a model fall within the reasonable 
ranges defined by simulations. These results assert NFW profiles as an effective model for the mass 
distributions of observed clusters. 

Subject headings: dark matter - clusters: individual (3C 220, A 370, CI 0016, CI 0024, CI 0054, 
CI 0939, CI 1409, CI 2244, MS 0451, MS 1137, MS 2137) - gravitational lensing 



1. INTRODUCTION 

While numerical simulations of dark-matter halos in 
the CDM model of cosmic structure formation invariably 
predict density profiles which are steeper than isother- 
mal outside and flatter inside a scale radius which is of 
order 20% the virial radius for cluster-sized halos (e.g. 
[N avarro ct al. 1997: Moore et al. 1998; Power et al. 2003; 
iNavarro et all2do3) , it is yet unclear whether real galaxy 
clusters h ave such density pro files. Galaxy rotation 
curves Csee lSofue fc Rubini2001lfor a review ) and strong- 
lensing constraints (e.g. Rusin & Ma '2001', ' Rusin et al.l 
12003; Trcu & KooDmans 2004; Kecton 2001) have shown 
that galaxies need to have at least approximately isother- 
mal density profiles which are however the result of bary- 
onic physics such as gas cooling and star formation. In 
galaxy clusters, baryonic effects should be substantially 
weaker, and thus their density profiles outside the inner- 
most cores should still reflect the typical CDM density 
profile found in numerical simulations. 

Gravitational lensing has been used in its strong and 
weak variants for constraining the density profiles of clus- 
ters. Weak lensing measures the gravitational tidal field 
caused by mass distributions, and thus allows density 
profiles to be directly inferred. While there is agreement 
among most studies of weak cluster lensing that clus- 
ter densi ty profiles arc compatible with the shape pro- 
posed bv iNavarro et al. (1996, 1997) (hereafter NFW), 
they are t ypically similarly well fit by isothermal pro- 
files (bm. IcWe et a1.l l200?l ICInwe Schueid^ l200Tt 
iSheldon et aLl 120011: lAthreva et al.ll2002() . The reason 
is that most of the weak-lensing signal comes from the 



cluster regions which surrounded the scale radius if the 
clusters had NFW density profiles, and there the NFW 
profile has an effective slope close to isothermal. 

Strong lensing can happen in the cores of sufficiently 
dense and asymmetric clusters and gives rise to highly 
distorted, arc-like images. There are now well over 60 
clusters known containing arcs with high length-to- width 
ratios. Mass models have been constructed for many of 
them, but mostly using axially-symmetric or elliptically 
distorted mass models with isothermal profiles. Several 
of these isothermal models turned out to be sp ectac- 

1 11991 IKneib et a,li ir996V 



ularly successful 
Constructed based on few large arcs, they were detailed 
and accurate enough to predic t counter-imag e s of a rclets 
found close to critical curves. iGavazzi et all |2003) find 
that t he core of MS 2137 seems to be closer to isothermal, 
while IKneib et alJ (2003) give an example for a cluster 
which is better fit by NFW than isothermal mass com- 
ponents. 

While models of strong lensing in clusters thus 
tend to favor density pro files steeper than expected 
from numerical s imula tions, iSand et all (|^004) followed 
iMiralda-Escudel l|1995fl in combining the location of ra- 
dial and tangential arcs with velocity-dispersion data on 
the central cluster galaxies and showed that cluster den- 
sity profiles should be substantially less cuspy in their 
cores than even the NFW profile. This conclusion hinges 
on the assumption of axial cluster symmetry and can 
be show n to break down for even mild ly elliptical mass 
models l)Bartelmann fc Meneghet!till2004fl . However, the 
situation is obviously puzzling, and it seems appropriate 
to ask whether samples of arc clusters can be success- 
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fully modeled with appropriately asymmetric NFW mass 
components. This entails two questions; first, can clus- 
ter mass models based on mass components with NFW 
density profiles be found which reproduce the observed 
arcs; and second, are the best-fitting model parameters 
within reasonable ranges defined by simulations? 

As our sample, we choose clusters which are known to 
have arcs and which have been imaged by the Hubble 
Space Telescope (HST). Our sample consists of f f clus- 
ters: 3C 220, A 370, CI 00f6, CI 0024, CI 0054, CI 0939, 
CI 1409, CI 2244, MS 0451, MS 1137, and MS 2137. We 
model each cluster with one or more elliptical NFW ha- 
los, each of which is completely defined by its scale con- 
vergence, scale radius, ellipticity, and position angle. In 
determining the values of these parameters that best re- 
produce the observed arcs, we constrain the mass distri- 
bution of the cluster. We find that all 11 clusters can be 
successfully modeled using elliptical NFW cluster mass 
profiles, and we tightly constrain the parameters defining 
each cluster's mass distribution. 

The rest of this paper is organized as follows: In § [3 
we describe our method for estimating the parameters 
describing a cluster's set of dark matter halos. In §01 we 
define our error estimation for the derived parameters. 
In § 01 we test our method of estimating lens parameters 
on a simulated lensing cluster of known properties. In 
§ O we outline our method of calculating the masses 
of our sample clusters. In § El we present results for 
each of the 11 clusters. Finally, in § 7, we discuss our 
main results and then summarize the implications of this 
work. Throughout this paper, we adopt a spatially flat 
cosmological model dominated by cold dark matter and 
a cosmological constant (firno = 0.3, Qao = 0.7, h = 0.7). 

2. LENS PARAMETER ESTIMATION 

Approximately 85% of the matter in galaxy clusters 
is dark. If at all, gas cooling plays a substantial role 
only in their innermost cores, where the gas density may 
be high enough for cooling times to fall below the Hub- 
ble time. It is yet unclear what influence gas physics 
may have on strong cluster lensing. While adiabatic gas 
seems to have little effect, efficient cooling and star for- 
mation may steepen the density profile very near the 
cluste r center and thus increase strong-lensing cross sec- 
tions l)Puchwein et al.i r2005'). In detail, any theoretical 
treatment of baryonic physics on strong cluster lensing 
depends on the numerical and artificial viscosity of the 
gas flow, the assumed star-formation efficiency, and the 
combination of a variety of feedback mechanisms. For 
simplicity, we shall here model each cluster in our sam- 
ple as a combination of purely dark matter halos. 

We model each dark matter halo with an asymmetric 
NFW profile. The spherical NFW density profile is 
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where ps is a characteristic density and Vs is the scale 
radius, which describes where the density profile turns 
over from p cx to p cx r~^. 

Following the common thin-lens approximation, the 
lens is approximated as a mass sheet perpendicular to 
the line-of-sight, and the scale convergence is defined as 
the ratio of surface mass densities, Hg = Ps^s/^crit, where 



Scrit is the critical surface mass density. 



(2) 



with the angular-diameter distances I?i,s,is from the ob- 
server to the lens, to the source, and from the lens to the 
source, respectively. 

Obviously, Kg is valid for a single source redshift only. 
In clusters showing arcs at multiple redshifts. Kg needs to 
be adapted in the following way. Assuming two source 
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> Zs"' for simplicity, we refer Kg 

to the lower source redshift zi^\ Fitting the lens model 
to the data, we adapt Kg by the factor 
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for the more distant sources, where D 
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Dis(2:i, zi''*), with zi as the lens redshift. 



(3) 

and 
Anal- 



ogous factors are applied for sources at additional red- 
shifts, if there are any. 

To elliptically deform the mass distribution, we alter 
the potential to have ellipsoidal rather than axial sym- 
metry. We accomplish this by defining the surface mass 
density E to have a radial dependence described by the 
projected elliptical radius. 



= [(r cos 61)2(1 - e) -f (rsin6l)V(l - e)] 
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rather than the circular radius r. We define the ellipticity 
e — 1 — b/a, where a and b are the major and minor axes 
respectively, and we define the position angle 9 in degrees 
counterclockwise from the +y axis. 

For a given cosmology and halo redshift, an ellipti- 
cal NFW halo depends on only four parameters: the 
scale convergence Kg, scale radius rg, ellipticity e, and 
position angle 0. Numerical simulations predict that the 
halo concentration, i.e. the ratio between the virial radius 
7^200 and the scale radius rg, is determined by the halo 
mass, albeit with considerable scatter (Navarro ct ^ 
T997; Bul lock et all Eimi lEke et al.1 120011: iDolag et all 
12004)) . This implies that the two parameters ps and rg 
characterizing a spherically symmetric NFW halo are not 
independent. However, with the aim of testing numerical 
results using strong cluster lensing, we do not adopt any 
correlation between these two halo parameters. 

We identify the arcs on the HST image of a cluster and 
define them by an array {xi,yi) of x and y positions of 
the image points constituting each arc. These points are 
arranged on a grid with spacing ct in a; and y. Generally, 
we take the spacing to be ct = 5 HST pixels in order 
to limit the number of arc points. This array of grid 
positions forms the data set we use to define a cluster's 
arcs. 

We use SExtractor l|Bertin fc Arnout£|ll996|) to define 
the position of each lens as its center of light on the 
HST image of the cluster. Then, for a given set of lens 
parameters (ks, Ts, e, 6) for each lens, we use the lensing 
equations to map the arc data back to the source plane. 
This yields a set of points describing the source. We next 
use the lensing equations again to map the source points 
back to the lens plane by finding all images of all source 
points. The result is a set of points that defines the image 
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of the source reproduced by the lens model. Our goal is 
to find the particular values of the lens parameters that 
yield predicted arcs that most closely match the arc data, 
the number of sources predicted by observations, and 
reasonably sized sources. Note that our approach does 
not require multiple arc-like images of a single source to 
be present and identified. 

For this purpose we define a function of the lens 
parameters to quantify the goodness of fit of our model 
to the data, so that the minimum of produces the best 
fit. Our x^ consists of three components. 

First, we define how well each data point is fit by the 
image points by finding the image point that lies closest 
to a given data point. We introduce a component Xi 
that depends on the distance between each data point 
and its closest image point. If we let the N arc data 
points be {xi,yi) and the closest predicted image point 
to a given data point be (ud^i, Vd^i), then this component 
of the x^ is 



to x^ from the source configuration as 



1 



N 
i=l 



(5) 



This Xi thus implicitly assumes that the reproduced im- 
age points are distributed in a Gaussian fashion centered 
on the data points, with a standard deviation of a. 

Second, we define how well each image point is pre- 
dicted by the data points by finding the data point that 
lies closest to a given image point. This is not redun- 
dant to the Xi calculation above because even if every 
data point has a nearby image point, there may still be 
distant, errant image points which are close to no data 
points. We thus introduce another x^ component xi that 
depends on the distance between each image point and 
its closest data point. If we let the M predicted image 
points be {uj,Vj) and the closest data point to a given 
image point be {xci.j,yci,j), then this component of the 
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assuming again a Gaussian distribution of data with re- 
spect to image points. 

Third, we require that arcs which observations indicate 
arc images of the same source indeed belong to a single 
source in our model. Some clusters host several fami- 
lies of arcs, each of which belongs to one unique source. 
We require that our model predicts both the number of 
sources suggested by observations of a cluster's arcs as 
well as the correct correlation between individual arcs 
and sources, as suggested by observations. In addition, 
we require that the predicted source be small and com- 
pact, as lensed sources are commonly observed to be. 

Define as the number of sources suggested by ob- 
servations to produce a given cluster's set of arcs. We ex- 
amine the zth source and its corresponding source points 
and image points predicted by our best-fit model to the 
cluster lens system. If we let the Pi predicted source 
points be (pij , qi,j)^ the mean pi position be pi, and the 
mean qi position be qi, then we define the contribution 




(7) 

where the source points are assumed to have a Gaus- 
sian distribution with standard deviation CTg. The choice 
of (Is is delicate as it controls the relative weight of the 
constraints in the lens plane, quantified by Xi 2i ^^'^ 
the source plane, quantified by xi- Large values of Cs 
may yield best fits with unreasonably large source con- 
figurations, while low values of (Tg may enforce very small 
sources at the expense of considerable deviations between 
image and data points. Thus, it may be necessary to 
try several fits with different choices of CTs to obtain both 
good agreement between image and data points and com- 
pact sources. 

We add these three components to yield our total x^, 

(8) 



2 2 
X =Xl 



xl 



xl 



By minimizing this quantity we maximize the overlap 
between data points and predicted image points, repro- 
duce the number of sources indicated by observations, 
and require that the sources be reasonably small in size. 
Each combination of (ks, Tg, e, 0) for each lens describes 
a different lens system that produces a different set of 
source and image points and hence a different x^- We 
use a downhill-simplex minimization routine ("amoeba" 
from lPress et al.lll992D to determine the combination of 
each lens's {ks, rg, e, 9) that minimizes x^- These best- 
fit values are the ones we use to define the cluster's lens 
system. 

3. ERROR ESTIMATION 

To estimate errors on the cluster dark matter halo pa- 
rameters deriv ed according to SEl we employ x^ statistics 
as outlined in iPress et alJ l)1992() . We must do so with 
caution. Formally, the x^ function is the log-likelihood 
assuming Gaussian distributions of model points relative 
to data points. We cannot be sure that this accurately 
describes our situation, in which we need to quantify 
the deviation between given data and reproduced image 
points. Assuming Gaussian likelihood factors, the two 
contributions Xi 2 quantify the likelihoods of reproducing 
the data points with the image points, and of finding im- 
age points exclusively near data points. Taken as another 
contribution to x^ , xi quantifies the log- likelihood of the 
sources being well modeled as Gaussians with widths CTg. 
Strictly speaking, our x^ is a figure-of-merit function as 
it quantifies the deviation of the model from the data in 
a well-defined sense. We shall, however, interpret it as a 
true x^ function, i.e. assuming that all deviations can be 
described by appropriate Gaussians. 

Recall that for a given lens system, the x^ function 
reaches a minimum for the best-fit parameter values. 
Since its gradient vanishes at the minimum, the x^ func- 
tion can to lowest order of a Taylor expansion be de- 
scribed as parabolic in the neighborhood of the mini- 
mum. By varying each of the parameters around this 
minimum, we thus expect to follow a parabolic section 
along the parameter axes through the x^ surface. 

For a cluster with iVi lenses, each modeled 
with an ellipsoidal NFW density profile, we 
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define a vector of best-fit parameters abf = 

(Ks,bf,l! '^s,bf,l, ebf,l, dhf,l, ■ ■ ■ , Ks,bf,Ni, fs,bf,Ni, ebf,Ni, ^'bt,N,)- 

We then calculate the 16 iVj^ components of the curvature 
matrix aki by noting that 

1 

2 



aki 



(9) 



where k and I vary from 1 to 4iVi. The partial deriva- 
tives can be calculated from the parabolic fits along all 
parameter axes to the local function. The covariance 
matrix is the inverse of the curvature matrix, C = a~^, 
and the 1-a error in each parameter is Suk = y/CJkk- We 
generally find these errors to be of order a few percent. 
An interesting approach for estimating errors on strong- 
lensing model parameters based on Monte-Car lo Markov 
Chains was recently proposed bv lBrewer fc Lewis (2Q05j) . 

4. COMPARISON OF PARAMETER ESTIMATION 
METHOD WITH SIMULATIONS 

To test our method of lens parameter estimation out- 
lined in §121 we use a numerically simulated galaxy cluster 
producing arcs. 

The cluster was kindly made available by Klaus Dolag. 
It was obtained by re-simulating at higher resolution 
a patch of a pre-existing large scale numerical simula- 
tion of the ACDM model with parameters ilmo — 0.3, 
^AQ = 0.7, Hq — TOkms""'^ Mpc~^ and normalization 
as = 0.9 . The "ZIC" techni que used is described in 
detail in iTormen et alJ l)1997j) . The cluster has red- 
shift z = 0.3 and a virial mass of Mvb = 2.29 x 
10^^ h^^ Mq. The particle mass in the re-simulation is 
Wpart 1-3 X 10^ h^^ Mq. The gravitational softening 
is set to 5h^^ kpc. 

We simulate lensing by this massive cluster using stan- 
dard ray-tracing techniques. First, the particles con- 
tained in a cube of 3ft~^Mpc comoving side length 
are selected. Then, to produce a two-dimensional den- 
sity field, their masses are projected along the line of 
sight, interpolating their positions onto a regular grid 
of 256 X 256 pixels using the Triangular Shaped Cloud 
method (Hocknev & Eastwood 1988). This surface den- 
sity map, shown in Figure ^ is used as the lens plane in 
the following lensing simulation. 

A bundle of 2048 x 2048 light rays is traced from the 
observer through the central quarter of the lens plane 
and their deflection due to the cluster mass distribution 
is c alculated as described in several earlier papers (see 
e.g. IMeneghetti erani2000l l200ll l2003b»al . The arrival 
positions of the light rays on the source plane, which we 
place at redshift Zg = 2, are used to reconstruct the lensed 
images of several sources distributed around the caus- 
tic curves so as to produce strong lensing features. The 
sources are modeled as ellipses, with random orientation 
and axial ratios randomly drawn with equal probability 
from [0.5,1]. Their equivalent diameter (the diameter 
of the circle enclosing the same area as the source) is 
= 1". Several arc configurations have been used to 
test our method, and some of these are illustrated in 
Figure 13 

We conduct a blind test of our parameter estimation 
method by applying it to the simulated cluster and arcs 
without knowledge of any of the cluster's physical prop- 
erties but its position. We permit knowledge of the clus- 
ter's position because when we model an HST cluster, we 




Fig. 1. — A 370" x 370" surface density map of the numerical 
cluster used to test our lens parameter estimation method. 
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Fig. 2. — The central region of Figure [Tl illustrating the critical 
lines for the numerical cluster as well as examples of arcs produced 
by this cluster lens. We used such systems of arcs to test our 
method of lens parameter estimation. 



determine the positions of its galaxy lenses on the HST 
image with SExtractor. 

We model the simulated cluster based on the cluster 
position and arc points, which is the same information 
we have when we model HST clusters. Applying the pa- 
rameter estimation method described in §|21 we estimate 
the scale convergence, scale radius, ellipticity, and po- 
sition angle of the simulated cluster. We then compare 
with the true values of these parameters in the simulated 
cluster mass distribution. 

In the simulation, the cluster is found to be well- 
described as an NFW ellipsoid with Kg — 0.54, rg — 
92.4", e = 0.18, and 9 = 99°. With our arc modeling 
method, we found the cluster's best-fit parameters to be 
Kg = 0.55, rg = 88.9", e = 0.17, and 9 = 96°, which are 
within 4% of the parameter values used to approximate 
the cluster. This is a convincing match, even more so 
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TABLE 1 

Best-fit parameters to lenses in the cluster, sample. 



Cluster 


Lens 




rs {h ^ kpc) 




e 


en 


CI 2244-02 




0.178 ± 0.003 


260 + 20 





113 + 0.005 


179 + 1 


Abell 370 


Gl 


0.164 ± O.OO?'' 


254 + 2 





28 + 0.01 


78 + 2 




G2 


0.165 + 0.0041" 


212 + 1 





073 + 0.004 


167 + 1 


3C 220.1 




0.178 + 0.002 


226 + 4 





265 + 0.005 


25 + 1 


MS 2137-23 




0.67 + 0.02 


64 + 2 





11 + 0.02 


95 + 3 


CI 140933+5226 




0.099 + 0.001 


241 + 2 





033 + 0.001 


79 + 2 


MS 0451.6-0305 




0.276 + 0.006<= 


262 + 8 





215 + 0.004 


156 + 1 


MS 1137+66 




0.256 + 0.002 


279 + 1 





143 + 0.005 


58 + 2 


CI 0054-27 


Gl 


0.047 + 0.001 


340 + 20 





175 + 0.003 


92 + 1 




G2 


0.100 + 0.001 


259 + 7 





14 + 0.01 


18.6 + 0.8 


CI 0016+1609 


DC 256 


0.102 + 0.003 


270 + 10 





123 + 0.002 


64.7 + 0.4 




DC 251 


0.087 + 0.001 


192 + 4 





121 + 0.005 


11.0 + 0.5 




DC 224 


0.219 + 0.004 


261 + 7 





1691 + 0.0008 


7.1 + 0.4 


CI 0939+4713 


Gl 


0.126 + 0.005 


136 + 6 





364 + 0.008 


2.65 + 0.03 




G2 


0.114 + 0.002 


190 + 10 





215 + 0.008 


73 + 1 




G3 


0.156 + 0.005 


170 + 1 





044 + 0.003 


87 + 3 


CI 0024+17 


#362 


0.158 + 0.003 


198 + 2 





059 + 0.002 


6.6 + 0.2 




#374 


0.170 + 0.002 


250 + 7 





153 + 0.006 


135 + 7 




#380 


0.116 + 0.003 


285 + 2 





0020 + 0.0001 


58.1 + 0.8 



''The value of Ks in the source plane of the giant arc AO. In the source plane of the arc pair 
B2/B3, Ks = 0.181 + 0.008 and in the source plane of the radial arc R, Ks = 0.24 + 0.01. 

'"The value of Ks in the source plane of the giant arc AO. In the source plane of the arc pair 
B2/B3, Ks = 0.183 + 0.004 and in the source plane of the radial arc R, Ks = 0.242 + 0.006. 

•^The value of Kg in the source plane of the upper arc, ARC2. In the source plane of the 
lower arc, ARCl, Ks = 0.55 + 0.01. 



TABLE 2 
Estimated cluster lens masses. 



Cluster 


Lens 


^arc 


rs M (< Ts) Reference 








(/i-ikpc) (/i-IMq) 



CI 2244-02 




2.237 


260 


9.33 X 10^^ 


1 


Abell 370 


Gl 


0.724/0.806/1.3 


254 


1.31 X lOi-* 


2 




G2 


0.724/0.806/1.3 


212 


9.21 X lO" 


2 


3C 220.1 




1.49 


226 


7.65 X 10" 


3 


MS 2137-23 




1.501 


64 


2.3 X 10" 


4 


CI 140933+5226 




2.811 


241 


3.78 X 10" 


5 


MS 0451.6-0305 




0.917/2.911 


262 


2.38 X lOi-* 


6 


CI 0939+4713 


Gl 


3.98 


136 


1.52 X 10" 


7 




G2 


3.98 


190 


2.68 X 10" 


7 




G3 


3.98 


170 


2.92 X 10" 


7 


CI 0024+17 


#362 


1.675 


198 


4.75 X 10" 


8 




#374 


1.675 


250 


8.17 X 10" 


8 




#380 


1.675 


285 


7.21 X 10" 


8 



References. — (1) 'Small et al." n997H); (2^ IBezecourt et alJ 11999^: C3 l 
E)ta et al. 12000); (4) Gavazzi (2005); (5) Lubin et al] 1200nr : f6) IBorvs et alJ 120041^ : 
f7) iTrager et al.. U997.') : (8) .Broadhurst et al. (.200011 . 



because the true cluster parameters are the results of a 
fit and carry errors themselves. Hence, we proceed with 
confidence in our lens parameter estimation routine to 
model our sample of HST clusters. 



5. CLUSTER MASS ESTIMATES 

If we know the lens parameters of a cluster and the 
source and lens redshifts, we can determine the mass of 
the lens within a given radius. We assume the lens is 
described by an NFW density profile, where the charac- 
teristic density is the density at the scale radius, given 
by ps = Ks'S'mt/fa- If wc define x = r/rg, then the mass 



contained within a (three-dimensional) radius R is 

j-B. 

M{< R) : 



:47rI]critKsr„ 



Hl + y) 



y 



(10) 



where y = R/rg. The critical surface-mass density is 
given in equation jSJ. 

For each lens in a cluster, we calculate the projected 
mass contained within the lens's scale radius. We note 
that for the clusters in our sample without published arc 
redshifts (MS 1137, CI 0054, and CI 0016), we cannot 
estimate the lens masses. 
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Fig. 3.— A 20" x 20" section of the F555W WFPC2 image of CI 2244-02 (left) and the same section of the best fit (right). The plot 
on the right illustrates the reproduced giant arc (gray points), the predicted source (black points), and the lens's critical curves. The gray 
image points are scaled such that darker grays denote brighter points in the image. 




Fig. 4. — As Figure^] but enlarged to 30" X 30" and using a model with lens parameters 3-a greater than the best-fit parameters (right). 



6. PARAMETERIZATIONS AND MASS ESTIMATES OF 
THE CLUSTERS 

Here we present the best fit to each cluster based on 
the observed arcs, as described in § |21 We will discuss 
each cluster individually. For a summary of the best-fit 
parameters for all sampled clusters, see Table ^ Om 
mass results are summarized in Table EJ which includes 
the reference for each cluster's arc redshifts. 

• CI 2244-02: The cluster CI 2244-02 has a red- 
shift z — 0.33 and host s a spectacu lar tangential 
arc (Small ct al. 1997bi). iLvnds fc P etrosian (198^ 
discovered this giant luminous arc, which is located 
near the cluster center. The arc is a partial Einstein 
ring, and we take the lens to be the large galaxy 
seen in Figure |21 near where the Einstein ring is 
centered. 

Although the errors we calculated for the best-fit 
parameters are as small as half a percent, they are 
indeed realistic. In Figure 01 we illustrate what the 



predicted source and images would be for a lens 
defined by parameters 3-tT greater than our best-fit 
values to CI 2244—02. The predicted giant arc is 
much larger than that observed, and is broken into 
two sections. Also, two small images are produced 
near the source that are not seen in observations. 

Changing the lens parameters by only 3-cr signif- 
icantly alters the lensed images, which indicates 
that the lens parameters must be tightly con- 
strained around our best-fit values. We use this 
example to justify the small errors we find for 
CI 2244—02, and extend the argument also to the 
small errors we find for clusters that follow. 

• Abell 370: Th e rich cluster Abell 370 is at red- 
shift z = 0.375 l|Abdelsalam et alJ 119981) and has 
a bimodal mass distribution with two cD galax- 
ies. These two galaxies mark the centers of two 
mass components in our model, and we identify 
the northern cD as Gl and the southern cD as G2. 
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By the classification of lBezecourt et aTl l)1999l) , the 
giant arc is AO, the nearby radial and two tangen- 
tial arcs are R, B2, and B3, and the upper tangen- 
tial arcs are Al and A2. These six arcs, as well as 
the cD galaxy lenses, are visible in Figure [S] 

We assume Al and A2 to be at the same redshift 
as the dominant arc AO. They are likely to form 
a double image of the same source. Arcs B2 and 
B3 share the same redshift and are images of a 
separate source; R is an image of yet another source 
l)Bezecourt et al.l 11999^. Accordingly, our best-fit 
model produces four independent sources for these 
six arcs. 

Our best fit also produces three additional images 
that we did not identify in our initial arc sam- 
ple. They correspond to fuzzy patches on the 
HST image which may reflect actual images. The 
exact number of arcs in the cluster is unclear; 
IBezecourt et al.l (|1999) identify as many as 81 ar- 
clets. Further study would offer more insight into 
the complex nature of Abell 370. 



• 3C 220.1: The cluster containing the radio galaxy 
3C 220.1 is at reds hift z = 0.62 and hosts a giant 
arc l|Ota et al.112000'1 . The lens responsible for pro- 
ducing this arc is 3C 220.1, which is identifiable as 
the largest source in the HST image on the left of 
Figure El 

• MS 2137-23: MS 2137-23 is a rich cluster at 
redshift z = 0.313 tha t contains a gi ant luminous 
arc and several arclets (jGavazz ill2002l) . The cluster 
is dominated by a single cD galaxy, which we define 
as the lens. A dopting the classification system of 
iGavazzil l|200;^ . the giant arc is AO, the radial arc 
near the cD galaxy is Al, the arc to the left of the 
cD galaxy in Figure [3 is A2, and the two arcs to 
the right of the cD galaxy are A4 and A5. The lens 
and the five arcs are visible in the HST image of 
MS 2137-23 shown in Figured 

IGavazzil l|2005D also fit an NFW model to 
MS 2137—23, and the cUipticity found in that pa- 
per is within the error bars of our best-fit cUip- 
ticity. However, our scale radius is ~ 30% smaller 
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Fig. 7.— As FigureOl showing a 45" X 45" section of the F702W WFPC2 image of MS 2137-23. 




Fig. 8.— As FigureOl showing a 10" X 10" section of the F702W WFPC2 imago of CI 140933+5226. 



than, and our scale convergence is twice as large as , 
the corresponding values found in iGavazz 
Despite these discrepancies, both our findings and 
those of ^Gavazzi (200^ agree that the concentra- 
tion of the cD galaxy in MS 2137—23 is high. 

Our lens model predicts an additional arc, near 
the center of the cD galaxy, which was also dis- 
cussed and p ossibly detected by IGavazzi et al.l 
iGavazzi (2002) suggests that the arcs in 
-23 can be categorized into two differ- 



(1200? 

MS 2137- 

ent systems, each with a unique source. The arcs 
AO, A2, and A4 are images of one source, and Al 
and A5 are images of a separate source. We use 
these correlations to constrain the lens parameters, 
and Figure[7|illustrates the two unique sources pre- 
dicted by our best-fit lens model. 



CI 140933+5226: iFischer et all l)1998j) discov- 
ered a quadruple lens in the z = 0.46 cluster 
CI 140933+5226, also known as 3C 295. The Ein- 
stein cross, as well as the lens galaxy at its center, 
is visible in Figure |H1 



As expected, the four images are due to a single 
source at the center of the Einstein cross. 

• MS 0451.6-0305: The cluster MS 0451.6-0305 
is located at z = 0.55 and hosts two tangential arcs 
|Borvs et al. 2004). The two arcs, as well as the cD 
galaxy lens, are apparent in the HST im age in Fig- 
ure El Following the classification of Borvs et alJ 
l|2004|) . the lower arc in the figure is ARCl and the 
upper arc is ARC2. 

iBorvs et al.l l|2004l) note that the two arcs in 
MS 0451.6-0305 are at different redshifts and thus 
are not images of a single source. These two sep- 
arate sources can be seen in our fit, shown in Fig- 
ure H 

• MS 1137 +66: The lensing c luster MS 1137-^66 at 
z = 0.783 (jClowe et al.lll998 D has the highest red- 
shift of the clusters in our sample and hosts several 
faint arcs. We model five arcs, and because there 
are no published redshifts for the arcs, for simplic- 
ity we assume that they are all at the same redshift. 
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The lens in our model is the cD galaxy at the center 
of the cluster. 

The lens defined by the parameters in Tabled ac- 
curately reproduces the five observed arcs, as seen 
in Figure ^| Each arc has its own independent 
source. 

We cannot estimate the mass of the cD galaxy lens 
because there are no published redshifts for the arcs 
in MS 1137+66. 

• CI 0054-27: The cluster CI 0054-27 is located 
at redshift z = 0.56 and has one lensed arc 
l|Smail et alJll997a|) . To reproduce the arc, we take 
the two lenses in the cluster to be the central galaxy 
(which we call Gl) and the upper left galaxy (which 
we call G2) in Figure [TTI 

These two lenses combined reproduce the observed 
arc accurately. A mass estimate of the two lenses 
in CI 0054—27 is impossible because there is no 
published redshift for the arc. 



• CI 0016+1609: The cluster CI 0016+1609 at red- 
shift z — 0.545, has one thin lensed arc l)Lavervl 
Il996j) . The three approximately coUinear elliptical 
galaxies seen in Figure E| define the center of the 
cluster and are the lenses we use to reproduce the 
observed arc. These three giant galaxies are, from 
top to bottom in Figure [HI DC 256, DG 251, and 
DG 224. 

The combination of these three lenses reproduces 
the single thin arc convincingly. Because there is 
no published redshift for the arc, we are unable to 
estimate the masses of the lenses in CI 0016+1609. 

• CI 0939+4713: The cluster CI 0939+4713 is at 
redshift z = 0.41 and has three radial arcs near 
its center (jSeitz et al.iri996.') . Three giant elliptical 
galaxies are visible in Figure El and these galaxies 
make up the cluster's core. We take these galaxies 
to be the gravitational lenses, and call the upper 
one in Figure 1121 Gl, the leftmost one G2, and the 
rightmost one G3. The three arcs are also visible 
in Figure [TSl 




Defined by their best-fit parameters given in Ta- 
ble n tfie three lenses are able to r eproduce the 
three arcs observed in CI 0939-1-4713. iTraeer et al.1 
l|1997D argue that two of the arcs in CI 0939-^4713 
are likely images of the same source, whereas the 
third arc is probably an image of a separate galaxy. 
Consistent with this expectation, our model pre- 
dicts two unique sources for the arc system. 

• CI 0024-1-17: The cluster CI 0024-^17 has redshift 
z = 0.395 and hosts five images of a single back- 
ground galaxy f Ota et ani2004j) . The three cen- 
tral elliptical galaxies in the cluster, seen roughly 
collinear in Figure 1141 act as the lenses for this 
system. From left to right in Figure m these 
galaxies are labeled #362, #374, and #380 in the 
iCzoske et all 1)200 ID catalog. 

As lenses, these three central galaxies can ac- 
curately reproduce the five observe d arcs in 
CI 0024+17, as shown in Figure [H lOta et all 
12004} note that all five arcs are images of a single 
background galaxy, and our results are consistent 



with that observation. 

7. CONCLUSIONS 

We have modeled 11 clusters hosting gravitational arcs 
and observed with HST as systems of dark matter halo 
lenses defined by elliptical NFW density profiles. Given 
its position on the HST image, each lens is completely 
defined by its scale convergence, scale radius, elliptic- 
ity, and position angle. We use a minimization routine 
to vary these parameters for each lens until the repro- 
duced images match the observed arcs in the cluster, and 
each observationally confirmed family of arcs belongs to a 
unique source. We also require that the predicted sources 
are compact. With this routine, we find the best-fit scale 
convergence, scale radius, ellipticity, and position angle 
for each lens in each cluster. 

Our main results can be summarized as follows: 

1. Each cluster in our sample is successfully modeled 
as a system of mass components with asymmetric 
NFW density profiles. The model produces images 
that correspond to observed arcs in the cluster, and 




the system of arc sources suggested by observations 
is reproduced. 

2. The best-fit parameters to each modeled cluster fall 
within the range of reasonable values set by simu- 
lations. Also, our estimates of the lens masses are 
reasonable for large galaxies. 

3. The accuracy of our minimization technique for 
finding the best-fit parameters to a system of NFW 
ellipsoids is verified by our fit to a simulated cluster 
with arcs. The cluster was simulated as an NFW 
ellipsoid, and we successfully reproduced the simu- 
lated cluster's parameters to within 4%. 

Clusters containing arcs have been commonly and suc- 
cessfully modeled using mass components defined with 
spherical or ellipsoidal isothermal density profiles. We 



show, however, that galaxy clusters can also be convinc- 
ingly modeled with NFW ellipsoids. The NFW profile is 
predicted for dark matter halos in numerical simulations 
of cosmological structure formation, and our result lends 
more credibility to NFW profiles as models of actual ob- 
served galaxy clusters. 
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